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We consider a two-dimensional lattice system with two sites in its unit cell. In such a system, 
the Bloch band spectrum can have some valley points, around which Dirac fermions appear as 
low-energy excitations. Each valley point has a valley spin ±1. In the system, there are two 
topological numbers counting vortices and merons in the Brillouin zone, respectively. These 
numbers are equivalent, and this fact leads to a sum rule that states that the total sum of the 
valley spins is absent even in a system without time-reversal and parity symmetries. We can see 
some similarity between the valley spin and chirality in the Nielsen-Ninomiya no-go theorem 
in odd-spatial dimensions. 
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The Bloch band spectrum in some lattice systems 
have valley points, around which Dirac fermions appear 
as low-energy excitations. The Dirac fermion has at- 
tracted much attention over the past few decades, since 
it is closely related to the quantized Hall effect, 1_7 ' and 
also plays important roles in the argument for topo- 
logical insulators. 8 ' The low-energy electronic features 
of graphene 9 ' and the zero-gap organic conductor a- 
(BEDT-TTF) 2 I 3 10,11 ' can be well described by gapless 
Dirac fermions. 

In this study, we consider a two-dimensional lattice 
system with two sites in its unit cell and discuss valley 
points. Each valley point has a valley spin ±1, which 
is well-defined as long as the intervalley mixing can be 
neglected, i.e., in the long-wavelength limit. We focus on 
the relation among the valley spin and two kinds of topo- 
logical numbers counting vortices 12, 13 ' and merons 14, 15 ' 
in the Brillouin zone, respectively. It has been shown 
generally that these numbers are equivalent on a two- 
torus. 16 ' Using long-wavelength formalism, we show that 
this equivalence leads to the fact that the total sum of 
the valley spins is absent. This sum rule is obvious when 
the system preserves time-reversal or parity symmetry, 
since a valley spin flips under these symmetry transfor- 
mations. In our argument, we do not introduce any as- 
sumptions on these symmetries. Namely, we can show 
the sum rule more generally. We also discuss some simi- 
larity to the Nielsen-Ninomiya no-go theorem for lattice 
fermions in odd-spatial dimensions. 17 ' We use the natu- 
ral unit h = c = 1. 

We start from a tight-binding model on a two- 
dimensional lattice system with two sites in the unit cell: 
H = J2im limcjcm, where cj (q) is the creation (annihi- 
lation) operator of an electron at the l-th site, and ti m is 
the hopping parameter between the l-th and m-th sites. 
In this system, we can introduce the sublattice spins A 
and B with the associated Pauli matrices <x, and the 
Bloch Hamiltonian can be written as ri k = e k cr + d k ■ cr, 
where fc is the two-dimensional crystal momentum, rj is 



the 2x2 unit matrix, and dk is a three-dimensional vec- 
tor in the sublattice spin space. The real spin is omitted, 
i.e., we consider a system where the real spin is not ac- 
tivated. The eigenvalue is = e k ± |dfe|. We note that 
the condition for crossing between the upper and lower 
energy bands is d\ k = dik = d 3 k = 0; and therefore, 
the most probable crossing is pointlike touching at some 
points in the two-dimensional Brillouin zone. More rig- 
orous arguments are given in the two-band case 18 ' and 
general case. 19 ' We use the following expression for the 
eigenstate: 

. (±K _ 1 f d 3k ± \d k \ \ m 
1 k '~n£\ d lk +id 2k j' W 

where is the normalization. We note the fact that 
we may choose the other expression \u^') cx (— d\k + 
id2k,d 3k T \dk\) T - 

In such a model, there can be some valley points in the 
energy spectrum. Around the z-th valley point fcj, the 
expansion of d k is given by d ki+Sk = d k% + 5k x d kx dk z + 
Sk y d ky d ki + 0(Sk 2 ), where 

d ki ■ d kj: dk z = d kt ■ d ky d kz = (2) 

is satisfied [here, we introduced the notations 
(<9fc x ,<9 fey ) = (d/dk x ,d/dk y ) and {d kx d ki , d ky d ki ) = 
(d klc d k , d ky d k )\ k = kl ]. In principle, we can systematically 
find out all of the valley points when we obtain an 
explicit form of d*, from the tight-binding Hamiltonian. 

We show a procedure for obtaining the "canonical 
form" of d k , which leads to the Dirac equation directly. 
The argument presented here basically relies on that 
given by Oshikawa. 7 ' We consider the rotation d 7 ^ — > 
J2i d l k R l k n , where I, m = 1, 2, 3 denote the coordinates in 
the sublattice spin space. We see that the Hamiltonian 
is transformed to H k — > UkHkUj., where Uk is a SU(2) 
matrix satisfying J2 m R 1 ^®™ 1 = U k a l U* k . Then, the state 

vector is transformed to 1"^ ) — > U k \uk). Using this ro- 
tation, we can align d ki along the e 3 axes. It seems that 
we may direct dk z parallel or antiparallel to e 3 ; however, 
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such an ambiguity can be removed as follows. First, we let 
d ki — > \dki\e3. We focus on the lower energy state (the 
discussion for the upper energy one is completely paral- 
lel), and (A) if \u^) -)• (0 1 cx.p[ita,n" 1 (d2kJd lki )}) T , it 
is consistent with the behavior of expression (1) we have 



Let us examine the sum rule 



chosen. However, (£>) if \u k 



(0,1) 



it is inconsis- 
tent. This behavior coincides with that of the other ex- 
pressions for \u k . ^). Thus, we flip d ki — > — \d ki \e^ [i.e., 

rotate 7r around the e 2 axes ] and then \u k . } ) -» (-1, 0) T , 
which is consistent with cq. (1). After performing an ap- 
propriate transformation, we introduce the mass parame- 
ter rrii = d^ ki , which is positive in case (.A) , and negative 
in case (B). Obviously, \rrn\ gives the band gap at fc,. 

Now, the vectors dk x d ki and dk y d ki are in the plane 
perpendicular to e 3 [see, eq. (2) ]. These two are not 
orthogonal, in general. Let us define the rank- 2 ten- 
sor A^l = 9^4, • dk„d ki , where fi, v = x,y. This is 
a symmetric tensor. This tensor would be regular at a 
generic valley point. Thus, we can find out the principal 
axes (k x i , k v > ) given by the rotation of the original axes 
(kx-iky) [parity transformation (k x ,k y ) — > (—k x ,k v ) is 
prohibited here, since it changes the valley spin defined 
below] and then A^) u , becomes diagonal, i.e., dk x ,d ki _L 
dk y idki- Moreover, we rotate the axes (k x ',k y >) so that 
dk x ,d kz || ei and dk y ,d kz || e 2 . We define v ifJ ,> = \d kfj ,d ki \ 
as the velocity. Then we arrive at the canonical form 
[hereafter, we omit the symbol / for the principal axes] 



d ki+ s k = {TiVi X 5k x ,Viy8k y ,mi) + 0(Sk 2 ), 



(3) 



where Tj = ±1. This is the valley spin that denotes the 
relative sign between the coefficients of 5k x and Sky. 9 ^ 11 ^ 
Here, we use the unitary transformation U ki = 03 to 
set the coefficient of Sk y positive. We should note that, 
in the above procedure, the parity transformations for 
k and the sublattice spin at a point on BZ should be 
prohibited, since we can change the sign of Tj arbitrarily 
using these transformations. 

The canonical form eq. (3) leads to the Dirac equa- 
tion, and the excitation near ki is described by the Dirac 
fcrmion and becomes dominant when the band gap given 
by 2 1 mi I is sufficiently small compared with the band- 
width. The Dirac cone arises 9-11 ) in the limit \m,i\ — > 0. 
We also note that the vector eq. (3) takes the configura- 
tion of a meron. 14 ' 15 ' 



In the limit of 5k — > 0, 



ki+Sk/ 



(mi > 0), 
(mi < 0), 



(4) 



where Q^ k = tan -1 (vi y 5k y ) / (vi X Sk x ). Namely, the lower- 
band state has a vorticity equal to Tj = ±1 around ki 
with mi > 0. 

We will see that the valley point specified by eq. (3) 
gives elemental contributions to topological numbers, 
which will be discussed later. For simplicity, we assume 
that there is no point which gives higher contributions to 
the topological numbers. Such a point can be recognized 
as an overlap of valley points. 



0, 



(5) 



where the summation is taken for all of the valley points 
in the band. This rule is obvious when the system pre- 
serves time-reversal and/or parity symmetry, since the 
valley spin is odd under time-reversal and parity trans- 
formations. Actually, this rule is satisfied in symmetric 
systems like graphene 9 ' and the zero-gap organic con- 
ductor a-(BEDT-TTF) 2 I 3 . 10 < n ) 

We show the sum rule in a more general manner. We 
do not introduce any assumptions on time-reversal and 
parity symmetries. We consider two topological numbers 
for the lower energy band: 



■Mr 



BZ 



BZ 



d*k. 

2TTi 

cPk 

477 

,(-)\ 



V fe x A k , 



(6) 



di. 



dk x d k x di 



d k ) , (7) 



where A k = (u k '\V k \u k '), d k = d k /\d k \, and J BZ d 2 k 
denotes the integral over the entire Brillouin zone. These 
two numbers characterize the topological structure of the 
lower band. 20 ' Generally, we can show that 16 ' 



Avor — A/rr 



(8) 



Let us estimate these numbers using the long- 
wavelength formations eqs. (3) and (4). JV vor is the Chcrn 
number that counts the total vorticity of the lower-band 
state, 12 ' 13 ' i.e., 



d 2 k 



BZ 



2m 



7 V fc x A k 



E 



dk ■ A k 



(9) 



around ki 



Ti8(mi 



0,±1,±2, 



where 0(x) is the step function. We can show it as follows: 
We may write A k = A k or + A k . The first part A™ 1 ' is 
defined to pick up the vortex singularity of the lower- 
band state eq. (4). Namely, 

AI 01 = iV k n9(mi) tan' 1 ^ ~ ^ | , (10) 

I . Vi x yK x Ki x ) J 

which gives the bottom line of eq. (9). 13 ' The remaining 
part Ak is, therefore, regular and does not contribute to 
A/^or, since the integral is defined on the entire Brillouin 
zone (two-torus). 13 ' Then, we obtain eq. (9). 

On the other hand, A/" m er is related to merons in the 
Brillouin zone. 14 ' 15 ' First, we vary all of the mass param- 
eters {m,i\i = 1,2...} to be infinitesimal without closing 
the gaps, i.e., without sign changes. This deformation 
does not change A/" mor because of its topological nature. 
From the configuration of a meron shown in eq. (3), we 
see that the integrand around fej is 



d k ■ [dk^dk x d k d k ) = 



(vfjk 2 



vlJkl + m 2 ) 3/2 : 



(11) 



which is valid for a small 8k. To estimate its contribution 
to A/rnor, we should introduce the appropriate momentum 
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cutoff Ask around ki, which is comparable to the band- 
width on the energy scale. On the other hand, the inte- 
grand eq. (11) is localized strongly around k i: since its ex- 
tension is characterized by an infinitesimally small \rrii\. 
Therefore, the contribution around ki can be obtained 
from the integration of eq. (11) without a momentum cut 
off, which gives a half-quantized number, sgn(mi)T;/2. 
This is the so-called "parity anomaly" . 1 ~ 7, 14, 15 ' Then, 
■Mner is given by the sum of contributions from all valley 
points, i.e., 

-d k ■ (d k:c dk x dk y dk) = ^2 2 S S n ( m T i- (12) 



BZ 



47T 



Eq. (12) has been verified numerically for arbitrary val- 
ues of mass parameters using an explicit form of dk 
on the entire Brillouin zone in a certain tight-binding 
model. 21 ) 

We can see from eqs. (8), (9), and (12) that the number 
of valley points should be even. We suppose that there 
are two valley points fe = fci , fe 2 . The extension to a sys- 
tem with 2N valley points (N = 2, 3, 4 • ••) is straightfor- 
ward. We emphasize that time-reversal and parity sym- 
metries are not required here: We do not introduce any 
restrictions on the locations of valley points or on the 
values of the mass parameters. Below, we show that eqs. 
(9) and (12) give the same result and become consistent 
with eq. (8) when (I) the sum rule (5) is satisfied, but do 
not when (II) the sum rule is not satisfied. 

We examine case (I) first. We can put t\ = +1 and 
T2 = — 1 without losing generality. From eq. (9), we ob- 
tain 



T\ + T"2 = Vfl\ , m 2 > 0, 

T\ = +1 mi > 0, m 2 < 0, 

r 2 = —1 mi < 0, m 2 > 0, 

mi, m 2 < 0. 



(13) 



On the other hand, from eq. (12), we obtain 



AT* 



mi 
\mi\ 



TO 2 

|m 2 | 



r r 2 



mi, m 2 > 0, 

+1 mi > 0, m 2 < 0, 

— 1 mi < 0, m 2 > 0, 

mi, m 2 < 0. 



(14) 



Namely, 7V VO r = A/" m or- In case (II), we immediately see 
that A/"vor 7^ Mncr, which is inconsistent with eq. (8). 

To sum up, we have shown that two topological num- 
bers shown by eqs. (6) and (7) estimated in the long- 
wavelength formalism give a result consistent with the 
general relation eq. (8), when eq. (5) is satisfied. 

Let us discuss the configurations of vortices and 
merons, which are schematically shown in Fig. 1. We 
note again that the valley spins arc fixed at t\ = +1 and 
r 2 = — 1. We note that a vortex with a vorticity equals to 
Ti is located at ki with m, > [see eq. (4)], on the other 
hand, a meron with a fractional charge sgn(m,:)T,;/2 is 
located at each ki. Therefore, when (a) mi,m 2 > 0, a 
vortex and a meron are located at fei, while an antivor- 
tex and an antimeron are located at fc 2 . In the case that 



(b) mi > 0, m 2 < 0, a vortex is present at ki but absent 
from fc 2 , while a meron is located at fci and fc 2 . In the 
case that (c) mi < 0, m 2 > 0, an antivortex is present at 
fc 2 but absent from fci, while an antimeron is located at 
fci and fc 2 . When (d) mi,m 2 < 0, there are no vortices 
in the entire Brillouin zone; however, a meron and an 
antimeron are located at fc 2 and fci, respectively. 





vortex 



anti-vortex /~\ meron V/ anti-meron 



Fig. 1. (Color onlone) Configurations of vortices and merons on 
the Brillouin zone (depicted by the two-torus). The valley spins 
are fixed at ti = +1 and T2 = —1. The signs of the masses are 
mi, 1712 > in (a), mi > 0, m2 < in (b), mi < 0,m2 > in 
(c), and mi,m2 < in (d). 



The phases (b) and (c) in Fig. 1 have nonzero topo- 
logical numbers and break parity and time-reversal sym- 
metries. 20 ) Such phases emerge if some interactions give 
mi and m 2 with opposite signs. Actually, this mecha- 
nism has been proposed by Haldane, 6 ' and its extension 
to the (time-reversal invariant) quantum spin Hall effect 
is given by Kane and Mele using a spin activating in- 
teraction. 22, 23 ) As long as the real spin is not activated, 
a system with parity and/or time- reversal symmetry is 
categorized into zero-topological number phase (a) or (d) 
in Fig. 1. We see that graphene 9 * 1 and the organic con- 
ductor 10,11 ' are the gapless limits of phase (a) or (d). 

The phase transition specified by the jump of the topo- 
logical numbers occurs when m^ changes its sign. 6,7,15 ' 
For instance, let us see a transition from phase (d) to 
phase (b) in Fig. 1, where the jumps of the topological 
numbers are AAf vor = AAA mcr = 1. We change mi con- 
tinuously from a negative value to a positive value. At 
the transition point mi = 0, a Dirac cone appears at 
fci in the energy spectrum, 6 ' and a vortex and a meron 
are created and an anti-meron is annihilated. The other 
Dirac-type spectrum at fc 2 with an opposite valley spin 
remains massive and hidden in the higher-energy part of 
the spectrum when its mass is comparable to the band- 
width. 6,24 ) 

Let us discuss some relations to the Nielsen-Ninomiya 
no-go theorem on the lattice fcrmion doubling in odd- 
spatial dimensions. 6, 17 > 18 - 24 ) When the system possess 
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parity symmetry, the sum rule eq. (5) becomes obvious 
as we mentioned earlier, and the number of valley points 
should be even. Besides, the energy gap at a point is for- 
bidden by the symmetry, since the mass term micr z comes 
from the symmetry-breaking staggered potential. (Note 
that, at a valley point, the parity-invariant mass term 
rrii<j z s z is allowed in graphene, 22 ' 23 ' but, we do not con- 
sider an interaction that activates the real spin) There- 
fore, we immediately see the fermion doubling. The role 
of the valley spin is similar to that of the chirality (the 
eigenvalue of 75 operator that can be defined in odd- 
dimensional space only) in the no-go theorem. A crucial 
difference is that the valley spin is still well-defined for a 
massive Dirac fermion as long as the intervalley scatter- 
ing can be neglected; on the other hand, chirality is not. 
Our discussion would be related to the argument given 
in rcf. 18 ' in which a two-dimensional analog of "chiral 
symmetry" is introduced artificially; however, the former 
appaers rather simpler. 

The situation becomes somewhat indefinite in parity- 
symmetry-breaking systems. As we have shown, the sum 
rule (5) is also satisfied in such systems. Owing to the 
sum rule (5), if we found a valley point, there should be 
another point with an opposite valley spin. We assume 
the following: (A) Band gaps at these paired points \m\\ 
and \m2 \ are degenerate, i.e., |mi = \ni2 \ = m. The pres- 
ence of time-reversal symmetry is a sufficient condition 
for this degeneracy. 25 - 1 (B) The gap amplitude m is suffi- 
ciently smaller than the bandwidth. Under these assump- 
tions, we find the statement: A massive Dirac fermion 
is always excited with its doubling partner that has an 
opposite valley spin, when the system has particle-hole 
symmetry (i.e., = 0) and the Fermi level lies in the 
gap. Assumptions (A) and (B) would be, at least, ap- 
proximately, satisfied when the symmetry-breaking per- 
turbations are small. 

The surface of a three-dimensional time-reversal in- 
variant topological insulator 8, 26 ' 27 ' provides an excep- 
tional case for the sum rule (5). In such a system, the 
topological #-term exists in the bulk region as the hall- 
mark of the Z2 topological order, and has a surface 
term that coincides with the Chcrn-Simons term with 
a half- quantized Hall conductivity. This fact indicates 
that a single Dirac cone without a hidden partner in 
the higher-energy region exists at a time-reversal invari- 
ant point in the Brillouin zone for the surface state, 
since such a cone gives a half-quantized conductivity 
owing to parity anomaly 1-7 ' 14 ' 15 ' and matches with the 
presence of the Chern-Simons term. Thus, the sum rule 
fails. The topological connection between the bulk and 
boundary regions, the so-called "bulk-boundary corre- 
spondence", 8 . 26 ' 27 ) causes this unique situation. This is 
somewhat similar to the fact that the doubling part- 
ner mentioned in the no-go theorem 17 ' with opposite 
chirality is absent in the one-dimensional edge state of 
the quantized Hall effect where the bulk-boundary cor- 
respondence is also at work. 28 ' 29 ' 

In summary, we have pointed out that the equivalence 
of two topological numbers [see eqs. (6), (7), and (8)] 
leads to the sum rule (5) for the valley spin defined in 



the long-wavelength formalism. The sum rule is obvious 
when the system preserves time- reversal and/or parity 
symmetry, since a valley spin is odd under these sym- 
metry transformations. In this study, the sum rule has 
been shown independently of the presence or absence of 
these symmetries. Basically, the valley spin at each val- 
ley point is determined by the detailed structure of the 
lattice tight-binding Hamiltonian. It seems interesting 
that a rigorous rule comes from a topological argument 
that is independent of the details of the Hamiltonian. 
We also emphasize that, to close the valley spin, we can 
see an analog of the fermion doubling theorem in odd- 
dimensional space 17 ' 18 ' in a rather simple manner. 
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